Geometric depolarization in patterns formed by backscattered light 
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We formulate a framework for the depolarization of linearly polarized backscattered light based 
on the concept of geometric phase, i.e Berry's phase. The predictions of this theory are applied 
to the patterns formed by backscattered light between crossed or parallel polarizers. This theory 
should be particularly adapted to the situation in which polarized light is scattered many times but 
predominantly in the forward direction. We apply these ideas to the patterns which we obtained 
experimentally with backscattered polarized light from a colloidal suspension. 
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The transport of light through human tissues is one 
of the most promising technique to detect in a noninva- 
sive way for instance breast cancer. For medical imaging 
applications, it is important to extract the information 
contained not only in the intensity but also in the polar- 
ization of backscattered light. This is not easy in general 
due to the complexity of vector-wave multiple scattering. 
In this paper we study a simple experiment, in which 
polarized light is backscattered from a diffuse medium. 
In these conditions, one observes between crossed po- 
larizers a fourfold symmetry pattern which was first in- 
terpreted qualitatively by Dogariu and Asakura 0. Re- 
cently more quantitative approaches have been developed 
for Mie scatterers using Mueller matrices A rather 
good agreement has been found between the experimen- 
tal shapes of the patterns and the theoretically predicted 
ones 

In this paper, we propose an alternate approach, which 
is very simple to implement because it is not based on 
a vector radiative transfer method as used generally in 
the literature. Instead our approach is based on the no- 
tion of geometric phase, which was introduced by Berry 
jsl] in his interpretation of the experiments showing op- 
tical activity in an helically wound optical fiber 6J. So 
far, the concept has been mostly applied to quantum me- 
chanics and to field theory but has not been used in the 
context of the transport of polarization in random media 
except in the recent ref. , which we follow and extend in 
this paper. Before presenting our framework, we discuss 
the cross-shaped patterns using the Stokes formalism to 
make contact with previous work Q Q . 

The Stokes parameters are the elements of the vector 
I = {I,Q,U, V) which is defined with respect to a plane of 
reference containing the direction of propagation. If the 
plane of reference is rotated through an angle 0, the new 
Stokes vector is L{(I>)I where is a rotation matrix 



We used the scattering matrix S corresponding to 
a distribution of spherical scatterers. The incident light 
is linearly polarized, is described by a Stokes vector Iq 
and is normal to the medium. Since the medium has 
cylindrical symmetry, the S matrix is independent of (/). 
After going through an analyzer (with Mueller matrix 
A), the outgoing Stokes vector is A ■ L{(t)) ■ S ■ L{(j)) ■ Iq. 
This means that the outgoing intensity in the crossed 
polarized (resp. copolarized) channel is 



I± = \ (2^11 



^33 - ^22) - \ {S22 + S33) COS 40, (1) 



/|| = ^ (2^11 - ^33 + ^22)-^12 COS 20+^ (^22 + ^33) COS 4 

(2) 

corresponding to Stokes parameters I = I± + In and 



Note that according to Eq. the crossed- 



polarized pattern has a four-fold symmetry, whereas ac- 
cording to Eq.[21 we see that the copolarized pattern con- 
tains in addition a two-fold symmetry due to the term 
proportional to 512 4]. In the particular case, which 
is satisfied in multiple light scattering Q, when S = 
Siidiag{l, C, C, D), with C ~ S22/ Sn and D — S44/S11, 
Eas. Ill2l take the simple form 



-Iq (1 - Ccos4(f)) 



i/o(l 



Ccos' 



(3) 



(4) 



corresponding to outgoing Stokes parameters I = Iq = 
Sii and Q = CIq cos icj). Note that a cross is expected 
now in both polarization channels and that C measures 
the contrast of these patterns. We have shown here that 
the patterns follow from general properties of symmetry 
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independently of the order of scattering or of the degree 
of coherence of the source. 

Let us now discuss the geometric depolarization of 
linearly polarized light. For Rayleigh scattering, the 
(linear) polarization vector after scattering E' is E' — 
k' X (E X k'), in terms of the polarization vector before 
scattering E and the scattered wavevector k'. This im- 
plies that E evolves by parallel transport in the limit 
of small scattering angles, and diffuse on the sphere of 
wavevector directions until the memory of the polariza- 
tion has been lost. Akkermans et al. has shown that 
this leads to a depolarization with a characteristic depo- 
larization length £p equal to ^/ln(10/7) ~ 2M Re- 
cently £p was measured using polarization resolved DWS 
[ifij i , which confirmed Akkerman's prediction for Rayleigh 
scatterers, and which gave ip ~ £* in the limit of forward- 
peaked scattering 17 — > 1 in agreement with Monte Carlo 
simulations Here we assume forward-peaked scat- 

tering because it applies to many biological tissues, and 
because in this case there is a clear analogy between light 
scattering paths and semi-flexible polymers . Recently, 
this analogy has been put on a solid basis, by realizing 
that the Fokker-Planck (FP) equation which describes 
semi-flexible polymers can be derived from the radia- 
tive transfer equation in this limit |l2|. In the following, 
we carry further the analogy by discussing the degrees 
of freedom of twist (polymer) analogous to polarization 
(light scattering). 

Let us consider a path of light, which we assume to 
be normally incident on a semi-infinite random medium. 
Following ref. Q, we express the polarization vector E 
in a basis of two vectors (n, b) normal to the tangent 
vector u (if the path is regular enough, the Frenet frame 
is a possible choice) as shown in figure 1: 

E{t)^ci{t)n{t)+C2{t)h{t), (5) 

where t is a parameter which goes from to s along 
the path. Let us call (j), the angle between E and n at 
i = 0, so that C2(0)/ci(0) = tan0. Since the polarization 
evolves by parallel transport c'l = tc2 and C2 — — tci, 
where t denotes the torsion on the trajectory, as found 
many years ago by Rytov fl3|. In the backscattcring 
geometry, n{t = s) = — n(< = 0) and h{t = s) = h{t = 
0), therefore we find that the polarization vector at the 
end of the path is 

E(i = s) = -cos(0-|-f7(s))n(s)-|-sin(0 + f^(s))b(s), (6) 

where n{s) is a geometrical phase, equal to the opposite 
of the integral of the torsion between t = and t — s 
modulo 47r 0. In the analogy between a path of hght 
and a semi-flexible polymer, the twist of the path is zero 
for light (it would be non-zero only in chiral medium), 
and the writhing angle is precisely 51. This writhe is a 
real value since the path is open, and that value is equal 
to the algebraic area of a random walk on a unit sphere. 



with the constrain that the path goes from the north pole 
to the south pole in the backscattcring geometry. From 
Eq. El we find that the output intensity after the light 
has gone through an analyzer crossed with respect to 
the direction of the incident polarization, is proportional 
to sin^(20 -t- Q). Because the medium is random, this 
intensity must be averaged with respect to all paths : 

I± (R) = J P'{s, R)ds < sm^{2(t) + fl{s)) >, (7) 

where P'{s,R) is the distribution of path length for a 
given distance to the incident beam R, and < .. > denotes 
the average over paths of length s. After expanding the 
r.h.s of Eq. [7| we obtain the form of Eq. |2|with Io{R) — 
J P'{s,R)ds and the contrast is 

C{R) = J P'i-s, R)ds < cos (217(,s)) > . (8) 

The factor cos(2i7) in Eq.|Hl means that the contrast re- 
sults from grouping pairs of paths of opposite geometrical 
phases, and the sum over s means that the phases of any 
other paths are uncorrelated. Interestingly a similar de- 
phasing occurs in the theory of magneto-conductance of 
Anderson insulators 0| . 

To evaluate the distributions of 11 for fixed s, P{s, fl) 
shown in fig. 2, we use a Monte Carlo algorithm originally 
developed for semi-flexible polymers. Random paths are 
generated with an exponential distribution of path length 
with a characteristic step equal to the elastic mean free 
path £. The incident photons are normal to the interface, 
but when exiting the medium all outgoing angles of emer- 
gent photons are accepted. The paths can be generated 
for an arbitrary ratio of the transport mean free path to 
the elastic mean free path £* /£. The geometric phase are 
calculated by closing the paths on the momentum sphere 
with a geodesic 0- Because of this closure, the distri- 
bution of n for short paths s ^ is peaked at zero as 
is also found for planar random walks (Levy's law). For 
long path s ^ £*, the distribution of fl widens, until 
the polarization is completely lost. In this regime, the 
distribution P(s,r2) is gaussian as required by the cen- 
tral limit theorem (the variance which was quadratic for 
s ^ £* becomes linear for s ^ ^* as seen in fig. 2), which 
implies that < cos2il(s) > is a decreasing exponential 
function of s. In Fig. 3, we show the corresponding curve 
for the contrast of the pattern calculated from Eq. |H1 to- 
gether with experimental points, which we obtained by 
averaging the Stokes parameter Q of an image along two 
perpendicular directions thereby suppressing a possible 
contribution in cos(2(/)) present in Eq. |2| In the experi- 
ment, a colloidal suspension of latex particles of negligi- 
ble absorption (diameter 0.5^m, wavelength A = 670nm) 
was used and the sample was about 8.8£* thick Q. The 
value of the anisotropy parameter g in the simulation was 
chosen to match the experimental value g ~ 0.82. In this 



3 



figure, we see that the contrast decreases exponentiaUy 
as function of the distance R with a characteristic dis- 
tance of the order of £p ~ £* both in the theory and in 
the experiments in agreement with refs. 
central region of the pattern, low order scattering is dom- 
inant as confirmed numerically. This could explain the 
discrepancy between experiments and simulations in this 
region, since our model only treats low order scattering 
events in an approximative way. 

To conclude, we have developed in this paper a 
simple theoretical framework for geometric depolar- 
ization, which we have applied to the experimental 
backscattering patterns. The mechanism of geometric 
depolarization is general provided that the scattering 
is peaked in the forward direction. We hope that our 
work will motivate further theoretical studies on the 
role of geometric phases in the transport properties of 
polarization in random media. 

We acknowledge many stimulating discussions with T. 
Maggs, M. Cloitre, F. Monti, and B. A. van Tiggelen. 
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List of Figures 

Fig. 1: Representation of a typical path in a semi- 
infinite random medium in backscattering. The Frenet 
frame consists of the tangent u, the normal n, and the 
binormal b vectors. R denotes the distance between end 
points, is the initial angle between the polarization 
vector E and the normal n, and 17 is the geometric phase. 

Fig. 2: Distribution of the geometric phase fl for 
different values of the path length s and in the inset 
variance of the distribution as function of s/£*. 

Fig. 3: Contrast as function of R: the crosses have 
been obtained from Monte-Carlo simulations using Eq.|H| 
and the squares are experimental values, obtained from 
an analysis of the Stokes parameter Q. 
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FIG. 2: Distribution of the geometric phase fl for different 
values of the path length s and in the inset variance of the 
distribution as function of s/£*. 



FIG. 1: Representation of a typical path in a semi-infinite 
random medium in backscattering. The Frenet frame consists 
of the tangent u, the normal n, and the binormal b vectors. 
R denotes the distance between end points, (j> is the initial 
angle between the polarization vector E and the normal n, 
and Q is the geometric phase. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

RJi* 



FIG. 3: Contrast as function of R: the crosses have been 
obtained from Monte-Carlo simulations using Eq. |H| and the 
squares are experimental values, obtained from an analysis of 
the Stokes parameter Q. 



